Irreversible magnetization under rotating fields and lock-in eff'ect on ErBa2Cu307_5 

single crystal with columnar defects 
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We have measured the irreversible magnetization (Mi) of an ErBa2Cu307_i single crystal with 
columnar defects (CD), using a technique based on sample rotation under a fixed magnetic field H. 
This method is valid for samples whose magnetization vector remains perpendicular to the sample 
surface over a wide angle range - which is the case for platelets and thin films - and presents several 
advantages over measurements of ML(-ff) loops at fixed angles. The resulting Mi(0) curves for 
several temperatures show a peak in the CD direction at high fields. At lower fields, a very well 
defined plateau indicative of the vortex lock-in to the CD develops. The H dependence of the lock-in 
angle ifL follows the H~^ theoretical prediction, while the temperature dependence is in agreement 
with entropic smearing effects corresponding to short range vortex-defects interactions. 



I. INTRODUCTION 

The study of the angular dependence of vortex pinning 
in high temperature superconductors (HTSC) with tilted 
columnar defects has revealed a richer variety of phenom- 
ena and pinning regimes than originally expected. At 
high temperatures and magnetic fields, the uniaxial na- 
ture of pinning by CD dominates the vortex responsecl. 
This is clearly seen, for instance, when isothermal mag- 
netization Iftops M(H) are measured for different field 
orientationaj. At fixed field modulus H, the irreversible 
magnetization Mi = ^AM (where AM is the width of 
the hysteresis, proportional to the persistent current den- 
sity J) exhibits a well defined maximum when H || CD. 
For other orientations "staircase pj/ortices" develop. In 
a previous study we have shownci that in YBa2Cu3 07 
(YBCO), and due to the simultaneous presence of CD, 
twin boundaries and crystallographic ab-planes, corre- 
lated pinning dominates over random pinning for all ori- 
entations, forming staircases of different configuration 
(i.e., with segments locked into different correlated struc- 
tures) depending on the field direction. 

An additional feature is the existence of a lock-in 
phase. When the angle between H and the CD is less 
than a lock- in angle (Pl{H,T), it is energetically con- 
venient for vortices to ignore ihe H orientation and to 
remain locked into the tracksB'Q. Since (pL scales as 1/H, 
in practice this effect is only visible at low fields. An 
experimaital manifestation of the lock-in regime is the 
existencea of a "plateau" in the irreversible magnetiza- 
tion. Mi (Q) « const, over a certain angular range. Here 
Q is the angle between the normal to the platelet crystal, 
n (which coincides with the crystallographic c-axis) and 
H, defined within the plane that contains the CD. 

At low fields, an additional effect must be taken into 
account. Due to both the anisotropic superconducting 
response of the HTSC and the sample geometry, the di- 



rection of the internal field B, that coincides with the di- 
rection of the vortices, differs from that of H. As the uni- 
axial pinning of the CD maximizes when B (rather than 
H) is aligned with the tracks, the maximum|-in Mi (O) 
occurs at an angle that progressively departsa from the 
orientation of the tracks, Od, as H decreases. 

The low field misorientation between B and H poses a 
serious experimental concern. All studies of the pinning 
properties of tilted CD that are based solely on measure- 
ments at H II CD, or on comparison of this orientation 
with a few others, give valid information at high fields, 
but are misleading at low fields; vortices are just not 
oriented in the right direction. To avoid this problem, 
a rather complete knowledge of the angular dependent 
response, Mi {H, 9), is required. 

In this work we present a procedure that allows us to 
obtain directly Mi (Q) by rotating the sample at fixed H 
and T. This method has the advantage that a fine grid 
can be easily obtained in the angular ranges of interest, 
thus permitting the exploration of the various regimes 
with significantly improved angular resolution. We apply 
this experimental procedure to investigate the pinning 
produced by tilted CD in an ErBa2Cu307 single crystal. 
We present a detailed analysis of the lock-in angle as a 
function of H and T. The width of the lock-in regime is 
shown to follow a 1/H dependence over a wide tempera- 
ture range, and from the temperature dependence of the 
slope of ipL vs 1/H we determine the entropic smearing 
function fiT/T*). 



II. EXPERIMENTAL DETAILS 

The sample used in this work is a rectangular 
ErBa2Cu307_5 single crystal platelet of dimensions 
0.44 X 0.33 X O.OlmTO^, grown by the self-flux rnethod in 
a commercial yttria-stabilized-zirconia crucibleB. After 



growth it was annealed under oxygen atmosphere for 7 
days at 450° C. This sample was then irradiated with 309 
MeV Am^^+ ions (whose penetration range in this mate- 
rial is ~ 15/i7TT,) at the TANDAR accelerator in Buenos 
Aires (Argentina), to introduce columnar defects at an 
angle Qd ~ 30° off the c axis. The rotation axis, which 
is perpendicular to the plane formed by the c-axis and the 
track's direction, is parallel to the largest crystal dimen- 
sion. The irradiation dose was equivalent to a matching 
field of i?$ = IT. After irradiation the sample presented 
a superconducting transition temperature of Tj, — 90.0K 
and transition width of AT < IK. 

Magnetization experiments were conducted on a com- 
mercial superconducting quantum interference device 
(SQUID) magnetometer (Quantum Design MPMS-5), 
equipped with two sets of detectors that allow to record 
both the longitudinal {Ml) and transverse {Mt) com- 
ponents of the magnetization vector M, with respect to 
the longitudinally applied field H. We have developed 
a sample rotation system (hardware and software) that 
solves the problems usually involved in the measurement 
of Mt, and thus allows us to study the response of the 
samples at arbitrary orientations. We have used that sys- 
tem in the past to measure magnetization loops ]V[(H) at 
different field orimtations in samples similar to the one 
investigated hercaEl In those experiments the sample is 
initially rotated to the desired 8, then zero- field-cooled 
(ZFC) from above Tc to the desired measuring temper- 
ature T. Both Ml{H) and Mt{H) are then measured, 
and the separation between the upper and lower branches 
of both loops {AMl{H) and AMriH)) is used to calcu- 
late the amplitude M, = i ^AMl{HY + AMt{HY and 
direction Qm = arctan{AMT{H) / /\Ml{H)) of the irre- 
versible magnetization vector Mj. 

The alternative rotating sample measurements pre- 
sented here are performed by setting up a desired ini- 
tial state {T,H,Qh) and then recording Ml{Q) and 
Mt(0) for fixed T and H. The sample is rotated a 
given angle step (typically 1° to 3°) and re-measured. 
Usually, the procedure is repeated until the crystal com- 
pletes 2 or 3 full turns. This provides us with redun- 
dant information that contributes to improve the qual- 
ity of the data. After careful subtraction of the sig- 
nal of the plastic sample holder (which has only lon- 
gitudinal component and is small, linear in H, almost 
temperature independent and, most importantly, angle 
independent) and of the reversible response, the irre- 
versible components MLi{H) and MTi{H) are used to 
determine M, = yj MLi{Hf + Mti{HY and Qm = 
arctan{MT^{H)/MUH))- 



III. ROTATING MEASUREMENTS 



A. Meissner response 

We began this study with an analysis of the Meissner 
response. To that end we ZFC the crystal, then applied 
a field H smaller than the lower critical field TJ^i (©) 
for all O, and subsequently performed the rotating mea- 
surements. Ideally, under those conditions there are no 
vortices in the crystal and the response depends neither 
on the material anisotropy nor on the pinning proper- 
ties, it is totally determined by the sample geometry. As 
was previously shownQ, Ml (O) and Mt (O) for a thin 
platelet should follow the dependencies 
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where v is the appropriate demagnetizing factor, which 
is essentially given by the thickness of the platelet {t) 
divided by its width {W) . These equations can be easily 
rewritten as 
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Equations y and [J indicate that the magnetization vec- 
tor M can be visualized as the sum of a fixed contri- 
bution Mo, anti-parallel to H, and a rotating contribu- 
tion M2© with a periodicity of 180°. This suggests that 
a convenient way to plot these data is on an Ml,Mt 
plane. In this presentation the Meissner response is ex- 
pected to lie on a circumference of radius M2e centered 
at {Ml,Mt) = {—Mo,0). One complete circumference 
is drawn by a rotation of 180°. An example of this pro- 
cedure (for T = 60K and H — 50Oe) is shown in Fig- 
ure 1. The crystal was rotated by two complete turns, 
thus there are four sets of data points covering 180° each, 
which are clearly separated in two groups. This is due 
to a small remnant magnetization Mr, which originates 
from the small residual field that is usually present during 
the ZFC.D 

The vector Mr has fixed modulus and its direc- 
tion rernains fixed with respect to the sample during 
rotation,Ll thus MfiL = Mfj cos (O -I- Qr) and M^t — 
Miism{Q + Qli), where Mn and e^; are constants. 



Since Mr (Q) has a one fold periodicity, it breaks the 
Meissner two fold periodicity and splits the experimen- 
tal data into two sets. Indeed, by fitting the data in 
Fig. 1 using a combination of the Meissner and remnant 
contributions we can easily determine and remove the 
remnant part and all points collapse on a single circum- 
ference (solid symbols). Fig. 1 is an extreme example 
of remnant influence, chosen to show that even in that 
case the Meissner response can be obtained. By carefully 
canceling the residual magnetic field in the ZFC proce- 
dure we can obtain a much smaller Mr such that both 
circumferences of raw data in Fig. 1 almost collapse on 
a single one. 

Equations and were used to fit several measure- 
ments for different temperatures and fields, and used to 
calculate the sample volume V Ki 1.45 x 10~^ cm^ and 
demagnetization factor i^ w 0.033. Both results are in 
very good agreement with the values directly determined 
from crystal dimensions (V ~ 1.46 x 10~'^cm'^; i^ ~ 0.03). 



B. critical state 

We now focus on the high field range, where the crystal 
is in the mixed state. Fig. 2(a) shows Mi(8) and Mt(0) 
for a rotation at 70K and 8kOe, where the angle inde- 
pendent background due to the holder has already been 
removed from M^. As the reversible magnetization of 
the superconductor [- ($o/327r2A2) \n{Hc2/H) - SG] is 
negligible compared to Mj, the response is dominated by 
vortex pinning. Curves in Fig. 2(a) exhibit a rich struc- 
ture, due to the combination of crystalline anisotropy, 
directional vortex pinning and geometrical effects. In or- 
der to extract useful information from them, we must 
first establish the relation between Mi and the screening 
current J flowing through the crystal. 

For simplicity, we will analyze the case of a thin infi- 
nite strip of aspect ratio v = t/W <C 1, that can rotate 
around its axis, which is perpendicular to H. Let's as- 
sume that the strip was originally ZFC at an angle Q 
and H was subsequently applied (the initial condition in 
Fig. 2(a)). If H is high enough we can considem that a 
current density of uniform modulus Jc{Q) flows over the 
whole volume.u This J is parallel to the strip axis and it 
reverses sign at the pJ^P^ that contains the axis and H. 

It has been shownETEil that, in this fully penetrated 
critical state and as long as vta-nQ ^ 1, the an- 
gle between Mj and the sample normal n is a ~ 
arctan (|i/^ tan8) <C Q- That is. Mi remains almost 
locked to n due to a purely geometrical effect. For 
the particular crystal of the present study, ajShould be 
smaller than 1° for 6 < 80°. Another result^ is that, 
although in principle the geometrical factor relating Mi 
with Jc{&) depends on 0, within that same angular range 



tated away from this initial state by a small angle SQ. 
The result will depend on the direction of rotation. If 
n approaches H (this corresponds to the angular ranges 
90° to 180° and 270° to 360° in Fig. 2(a)), the normal 
component H± will increase, thus inducing screening cur- 
rents at the edges of the crystal in the same direction 
as those already flowing. Vortices will then displace to 
satisfy the condition J < Jc{Q + SQ) everywhere. If 
Jc (6 -I- 5Q) < Jc(0) the new distribution will be anal- 
ogous to the initial one, with J — JciQ + SQ) every- 
where and the boundary of current reversal rotated by 
an angle SQ in order to remain parallel to H. On the 
contrary, if Jc(6 -I- SQ) > Jc(0), the new field profile 
will propagate all the way to the center of the sample 
only if SH± — H sm{Q)60, is larger than the maximum 
possible additional screening ~ t [Jc{Q + 5G>) — Jc(0)]. 
The condition for the "full penetration of the rotational 
perturbation" is thus 
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the variations are given by the factor (l — ^v^ ta.n^ Q) 
and thus are negligible. 

We now discuss what happens when the strip is ro- 



ll the inequality (R) is satisfied, the vortex system will 
evolve under rotations maintaining a fully penetrated 
critical state with uniform J. In other words, the state 
at any Q will be the same that would have formed by in- 
creasing H after ZFC at that orientation. Then, as long 
as i^tanO <C 1, the condition that Mi is almost paral- 
lel to n is preserved. We have experimentally confirmed 
this fact: a < 1° for all measurements conducted in this 
work, except in a very narrow angulatj-ange around the 
ab-planes, where a flip in Mi occurs.113 Thus, from now 
on we will plot all the results as a function of Q. Fur- 
thermore, we can obtain Jc(6) by simply multiplying Afj 
by the angle independent factor that corresponds to the 
relation valid for H || n. If eq. (1^) is not satisfied, J will 
be subcritical in part of the sample and this relation is 
no longer valid. 

If the crystal is rotated in such a way that n moves 
away from H (so H± decreases), the new screening cur- 
rents induced at the edges of the crystal will oppose to 
those already flowing. As the rotation progresses the 
boundary between the old and new J directions will move 
inwards, until eventually the new critical state propa- 
gates to the whole sample. From that point the situation 
will again be analogous to that already discussed, except 
that Mi will be paramagnetic instead of diamagnetic. 

A rotation at |£xed H is to some extent analogous to 
a hysteresis loopQ. Rotating n towards H increases H±, 
which is roughly equivalent to increasing _ff at O = 0°, 
moving along the lower (diamagnetic) branch of the loop. 
Decreasing H± (either by rotating n away from H or by 
crossing the H || c condition), is equivalent to revers- 
ing the field sweep, thus producing a switch to the other 
branch of the loop. This is a useful analogy for the anal- 
ysis of the rotations, although it should not be pushed 
too far. 

A basic difference is that a rotation also produces 
a variation in the parallel field component, SHu = 



H cos(8)(50. This generates screening currents flowing in 
opposite directions on the upper and lower surfaces ^ttlp^ 
strip, which produce a tilting force on the vortices.ll2rO 
If the perturbation propagates all the way to the cen- 
tral plane, the result is a rotation of the vortex direction 
following H, the situation that we have implicitly as- 
sumed above. However, if pinning were strong enough 
it could preclude the propagation of the tilt beyond a 
certain depth, thus generating a critical state along the 
crystal thickness, with a central segrae|nt of the vortices 
remaining in the original direction.EHrtS If this effect were 
significant, as the rotation proceeded the orientation of 
the vortices would lag behind the field direction. In an 
extreme case, vortices deep inside the sample would ro- 
tate rigidlifr-with it, a situation that has indeed been 
observedll3 ell. As we will show below, in the present 
case we have clear experimental evidence that the mis- 
orientation between the vortex direction and H due to 
this lag effect is negligible, so all this complication can 
be ignored. 

We now analyze the curves shown in Fig. 2. The mea- 
surement starts at 8 '^ 30° (point A) with n rotating 
away from H. Thus, J initially undergoes a flip until the 
reversed fully penetrated critical state is formed (point 
B). From here the evolution of the system turns indepen- 
dent of the initial conditions and becomes two fold peri- 
odic. From point C (6 := 90°) to point E (6 = 180°) the 
system evolves in a fully penetrated critical state (in the 
hysteresis loop analogy, this is equivalent to increasing 
the field from zero to H) . Clearly visible within this an- 
gular range is the peak in both Ml and Mt at 6 ~ 150° 
(point D), that corresponds to the direction of the CD. 
At point E, Mt is null as expected by symmetry, while 
Ml begins a quick flip due to the reversal of the screening 
currents as iJ_L reaches a maximum at H || n and then 
starts to decrease. The end of this flip at point G indi- 
cates that the critical state is completely reversed. From 
G to C (8 = 270°) the evolution is analogous to a field 
decreasing portion of a loop. 

Note that between E and G there is one unique angle 
(point F) where both Ml and Mt are null. This con- 
dition is equivalent to the unique H value in the switch 
from the lower to the upper branch of a M{H) loop where 
Mi = 0. The fact that the condition Ml = occurs at 
the same angle where Mt = confirms that the back- 
ground signal has been correctly subtracted, and we have 
systematically made use of this checking procedure. 

In Fig. 2(a) the direction of rotation is such that 
the conditions H _L c; H || CD and H || c proceed in 
that order. We define this as a clockwise (CW) rotation. 
In contrast, in a counter-clockwise (CCW) rotation the 
alignment occurs when n is moving away from H. The 
consequences of this difference are described below. 

In Figure 2(b) the same CW data of Fig. 2(a) is shown 
in an Ml vs Mt polar graph (full symbols) , together with 
the CCW rotation under the same conditions (open sym- 
bols). In both cases the initial behavior until the critical 
state is fully developed (portion A to B in the CW and P 



to Q in the CCW) and the subsequent 180°-periodic evo- 
lution in the critical state (covering approximately two 
periods of 180°) are clearly distinguished. Another fea- 
ture that is apparent in this representation is that the 
magnetization vector passes through the origin (Mi = 0) 
and reaches the opposite quadrant each time that (i) a 
rotation starts moving n away from H; or (ii) the n || H 
condition is crossed. 

Although the CW and CCW curves in Fig. 2(b) are 
similar (rotated in 180° with respect to each other) they 
also exhibit some differences. The most obvious one is 
that the peak at the CD direction (dotted line) is promi- 
nently seen in the CW rotation (point D), while in the 
CCW rotation it is partially suppressed by the flip of Mi. 
The flip starts at H || n, and ends at the angle Qp where 
the fully reversed critical state is achieved. Making use of 
the loop analogy, this requires a field decrease of ~ 2H*, 
where H*{H,T) is the well known full penetration field, 
then 



2H* = iJ(l-cos8F) 



(8) 



This analysis indicates that there is a blind range in 
the rotation measurements, extending up to an angle Qp 
from n, where the critical state is not fully developed and 
thus J cannot be extracted. Depending on the direction 
of rotation, this blind range occurs either in the same 
quadrant of the CD (case CCW) or in the opposite (case 
CW). As Qp decreases with H, in CCW rotations the 
peak due to the CD is totally hidden at low fields but 
can be fully measured at high enough H . 

The values of Op are easily obtained from Figure 3(a), 
where Mi is plotted as a function of 8 for the same 
two sets of data (CW and CCW) of Fig. 2(b). We ob- 
serve here that the agreement between the CW and CCW 
data is excellent, thus they can complement each other 
to eliminate the blind region at low angles. Estimating 
Sp - 25° for the CW rotation and 8f -- 30° for the 
CCW case, and using eq. (||) we obtain H* -- 370Oe 
and ~ 540Oe respectively. We can check the consistency 
of these estimates in two ways. First, we know that in a 
thin sample H* ^ Jt. Combining with the critical state 
relation J ~ 60Mi/W (valid for a square platelet) we 
have H* ~ 60M,t/W ^ l.SMi. From the figure we have 
MiiOp) - 200G for the CW and - 3400 for the CCW, 
so we get H* ~ 360Oe and ~ 610Oe respectively, in very 
good agreement with the above estimates. On the other 
hand, we can compare the values of H* obtained from 
eq. (^) with those directly measured in hysteresis loops 
at the appropriate angles. We have done so for several 
temperatures and fields, and we have systematically ob- 
tained very good consistency. 

Fig. 3(a) confirms that the condition (R) is satis- 
fied in this measurement. In fact, the largest slope 
dMi/dQ ^ IkG/rad, that occurs at 8 ~ 33°, implies 
that tdJc/dQ ~ 1.8kG/rad, which is indeed smaller than 
H sm{Q) ~ AAkG. This condition is also fulfilled in all 
the cases discussed in the next section. 



In order to compare the data measured by sample ro- 
tations with those resulting from traditional loop mea- 
surements, in Fig. 3(a) we also included Mi values at 
several obtained in the latter way at the same T and 
H (large open diamonds). The agreement is very good 
over the full range of angles, except that the loop values 
tend to be somewhat smaller. This is a feature observed 
for all measured fields, and can be explained by the fact 
that a rotation step is a process that takes only a couple 
of seconds, while a field increase and stabilization typi- 
cally requires more than 1 minute in our magnetometer, 
during which the Mi is already relaxing. Indeed, by per- 
forming short relaxation measurements we have verified 
that the rotation data approaches the loop data after 1-2 
minutes. This results in another advantage of the rota- 
tions over the loops: measurements are made closer to 
the true initial critical state. 

Finally, the coincidence of the CW and CCW rotations 
and the loops, particularly in the region of the peak due 
to the CD, rules out the possibility that vortices lag sig- 
nificantly behind the direction of H in our rotating sam- 
ple experiments. In summary, the information obtained 
from our rotation measurements is essentially the same 
as that provided by hysteresis loops, with several ad- 
vantages including the possibility to acquire significantly 
more data points for each field. This feature permits a 
more detailed analysis of the peak associated with the 
uniaxial pinning of the CD, as will be shown in the next 
section. 



IV. DETERMINATION OF THE LOCK-IN 
ANGLE 

A complete set of rotations at 7QK for different ap- 
plied fields is shown in Figure 3(b). At high fields (above 
~ 6kOe) a well-defined peak at the CD direction is ob- 
served. At lower fields (IfcOe < H < bkOe) the peak 
progressively broadens and transforms into a plateau (a 
certain angular range where Mi{Q) '^ const.), while ii, 
shifts towards the c-axis. We had previously reportedcl 
all these features in YBa2Cu307 crystals. 

The plateau represents the angular range of applied 
field over which it is energetically convenient for the vor- 
tices to remain locked into the columnar defects, thus its 
angular width is twice the lock-in angle (pL- Below this 
angle, the vortices are subject to an invariant (and ma 
imum) pinning force. According to theoretical models 
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where ei is the vortex line tension and Sr (T) is the effec- 
tive pinning energy. 

Equation (|9|) predicts that (p^ should be inversely pro- 
portional to H. The improved resolution of the rotation 
measurements, that permits a much better determination 



of the width of the plateau, allows us to test this depen- 
dence. To that end, we have measured several other sets 
of data similar to Figure 3(b), for a wide range of tem- 
peratures {35K to 85K). A few examples of the observed 
plateaus are shown in Figure 4. 

We then extracted the plateau width for every mea- 
surement which displayed such a feature. This proce- 
dure was done very carefully, including an over-zealous 
estimate of the errors involved. The results for all mea- 
surable (fL are plotted as a function of H~^ in Figure 
5. This figure clearly demonstrates the H^^ dependence 
of (/?!,, as evidenced by the solid lines which are the best 
linear fits to the data points for each temperature. 

According to eq. (||) , the data in fig. 5 should extrapo- 
late to the origin, what is clearly not the case. For all the 
temperatures where reliable extrapolations can be made 
(35-?^ to 80K) the linear fits systematically give a posi- 
tive value of ifiL ~ 1.5° to 3° at H~^ = 0, which is above 
the experimental error. There are at least two reasons 
for this discrepancy. In the first place, we experimen- 
tally determine ifL from the intersection of straight lines 
extrapolated from the plateau and the slopes at both 
sides of it (see fig. 4). Due to the rounded ends of the 
plateau, this definition tends to overestimate ipL- Sec- 
ond, the natural splay of the tracks will tend to wash 
away the expected cusp-like behavior at high fields, thus 
also contributing to the overestimate of Lp^. It is clear, 
on the other hand, that the influence of the splay is not 
too dramatic, as we indeed observe a rather sharp peak 
at high fi#Ws as seen in fig. 3(b). TRIM calculations 
show thatE3, in our irradiation conditions, the median 
radian angle of splay slowly increases from zero at the 
entry surface of the crystal to ~ 3° at a depth of %^m, 
and then grows faster to ^ 6° at the exit surface. 

We now want to analyze whether eq. (p|) provides a 
satisfactory description of the temperature dependence 
of the lock-in effect. As this expression does not ac- 
count for the nonzero extrapolation of Lp^ discussed in 
the previous paragraph, it would be incorrect to force 
a fit through the origin to determine the prefactor of 
H^^. Instead, it is appropriate to identify such prefac- 
tor with the slopes a{T) = d(fL/d{H^^) of the linear 
fits. Indeed, the splay of the CD is a geometrical fea- 
ture independent of H and hence it should only add a 
constant width to the plateau, without changing its field 
dependence. To a first approximation, the rounded edges 
of the plateau will also introduce an additive constant, 
without significantly affecting the slope. Figure 6 shows 
the temperature dependence of a{T) (solid symbols). As 
expected, Oi{T) decreases with increasing T, reflecting 
the fact that the lock-in angle at fixed H decreases with 
T due to the reduction of both the line tension and the 
pinning energy. For a quantitative analysis it is neces- 
sary to know the expressions for ei and er{T). In our 
experiments the appropriate line tension is that corre- 
sponding to in-plane deformations (see pages 1163-1164 
in Ref.a), ei = (£^eo/e(0)) Ihk, where Eq = ($o/47rA)^, 



the penetration depth A corresponds to H He, the mass 
anisotropy e < 1 and e^(6) = cos^(e) + £^ sin^(e). The 
temperature dependence of the superconducting param- 
eters appears in pgrf, through A(T). On the other hand, 
Sr{T) is given b; 
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where r sa SOA is the radius of the tracks, ^ is the su- 
perconducting coherence length, and the dimensionless 
efficiency factor 77 < 1 accounts for the experimental fact 
that the pinning produced by the CD is sniaUer than the 
ideal.E3 Besides the intrinsic temperature dependence of 
the superconducting parameters, this expression contains 
an additional temperature dependent factor /(x), known 
as the entropic smearing function, which accounts for the 
thermal fluctuations of the flux lines. Here x = T/T^p, 
where Tdp is a characteristic field-independent depinning 
temperature. Combining all these elements, at the track's 
direction O = Qcd we obtain 
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In the original work of Nelson and Vinokum, where 
only a short range pinning potential was considered, the 
entropic function for a: > 1 was approximately given by 
fsr{x) ~ x^exp(— 2x^).pHowever, according to a further 
refinement of the modelcl, where the long range nature of 
the pinning potential was taken into account, this func- 
tion (for X > 1) takes the form fir{x) ~ exp(— x). 

We can now fit the experimentally determined a{T) 
using eq. (|ll|). To that end we use the long range result 
fir{x) and fix the reasonably well known superconduct- 
ing parameters of the material e « 1/5; In k « 4 and ^ = 
15A/\/l — t (where t = T/Tc). We also assume the usual 
two-fiuid temperature dependence A(T) = Al/2\/1 — t*, 
where A^ is the zero-temperature London penetration 
depth. The free parameters are then Tdp and the combi- 
nation Ai/?7^/^. The best fit, shown in figure 6 as a solid 
line, yields Ai/r/i/" = 360A and Tdp = SOK. 

Based on the results of figs. 5 and 6, there are a num- 
ber of considerations that can be made at this point. The 
first one is that the Bose-glass scenario contained in eqs. 
(^ to (^l|) provides a quite satisfactory description of 
the lock-in effect over the whole range of temperature 
and field of our study. In addition, the obtained Tdp is 
smaller, but still reasonably similar to the value ^ AlK 
that we had previously found for several YBCO crys 
using a completely different experimental method.E^ 
This low Tdp (well below the initial theoretical expec- 
tations) indicates that the efficiency factor r] is rather 
small, what is also consistent with the less than opti- 
mum Jc observed here and in several previous studies. 
For low matching-fields as that used in the present work, 
it was estimateclE3 that 77 ^ 0.2 — 0.25. 



The exact value of r] has little influence in our estimate 
of Ai, as it only appears as rj^^^. For r] = 0.2 and rj — 1 
we get Xl = 250 A and 360A respectively, a factor of 4 to 
5 smaller than the accepted value Xl ^ 1400A. Zhukov 
et alxB had reported a similar discrepancy when study- 
ing the lock-in effect by both CD and twin boundaries in 
YBCO. In a previous study in YBCO crystals with CD, 
we had also found that the misalignment between B and 
H at low fields (due to anisotropy effects) was well de- 
scribed using a A^ significantly smaller than the accepted 
valua3. Thus, this numerical discrepancy appears to be a 
common result associated to the study of angular depen- 
dencies in YBCO-type superconductors with correlated 
disorder. 

Finally, it is relevant to note that eq. (y) was de- 
rived for the single vortex pinning regime, which occiixa 
below a temperature dependent accommodation fieldQa 
B*{T) < _B$, while a large fraction of the data shown in 
fig. 5 lies above this line, in the collective pinning regime. 
Unfortunately, to our knowledge there is no available ex- 
pression for ipi^ {H, T) in the collective regime. Blatter et 
al.Q only argued that collective effects should result in a 
reduction of the lock-in angle. The experimental fact is 
that eq. (||) satisfactorily describes both the temperature 
and field dependence of (pL- This suggests that, at least 
to a first approximation, collective effects in the range of 
our measurements simply result in a different prefactor 
in eq. (H) . Clearly, lock- in effects in the collective regime 
deserve further theoretical study. 



V. CONCLUSIONS 

We have measured the irreversible magnetization (Mi) 
of an ErBa2Cu307_5 single crystal with columnar defects 
(CD), using an alternative technique based on sample ro- 
tation under a fixed magnetic field. The resulting Mi(0) 
curves for several temperatures agreed very well with in- 
dependent hysteresis loop experiments, showing a peak 
in the CD direction at higher fields, while a very well 
defined plateau due to the lock-in of the vortices into the 
CD was observed at lower fields. The lock-in angle sat- 
isfactorily follows the field and temperature dependence 
predicted by the Bose-glass scenario. 
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Figure 1. Mt versus Ml polar graph in the Meissner 
phase for T — 60K and H — 50Oe. Open symbols: 
raw data. Solid symbols: pure Meissner response (the 
remnant contribution was removed). 

Figure 2. (a) Components of the magnetization vector, 
Mt and Ml, as a function of the angle for T = 7QK 
and H — 8kOe. The direction of rotation is such that 
the conditions H _L c; H || CD and H || c proceed in 
that order, (b) Mr vs Ml polar graph of the same CW 
rotation of (a) (full symbols) together with the CCW 
rotation. 

Figure 3. Irreversible magnetization M^ as a function 
of e at T = 70K for (a) H = 8fcOe together with data 
obtained from hysteresis loop measurements (b) several 
fields. 

Figure 4. Irreversible magnetization Mi as a function 
of O in the region of the plateau at T = 50K and 70K 
for several fields. 

Figure 5. Lock- in angle ifL versus l/H for several tem- 
peratures. The straight lines are fits according to equa- 
tion (|). 

Figure 6. Temperature dependence of the slopes 
ct{T) = dLpL/d{H~^) of the linear fits of fig. 5 (full sym- 
bols). The solid line is a fit to eq. ([ll]). 
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